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SUMMARY 


/l* 

In this work, numerical methods are developed for 
control of three-dimensional adaptive test sections . The 
physical properties of the design problem occurring in the 
external field computation are analyzed , and a design 
procedure suited for solution of the problem is worked out. 

To do this, the desired wall shape is determined by 
stepwise modification of an initial contour. The necessary 
changes in geometry are determined with the aid of a panel 
procedure or, with incident flow near the sonic range, with 
a TSP procedure. 

The poor convergence properties of the classical 
iterative adaptation process in the 3D case are demonstrated 
by means of a numerical simulation and a theoretical study. 

At the same time it is shown that by means of a simple 
replacement of the constant control factor by a matrix which 
considers the incident Mach number and test section 
geometry, computation of the adapted wall contour can be 
determined practically in one step, i. e., on the basis of 
a single wind tunnel test. The capability of this one-step 
method is demonstrated by means of measurements from the 
adaptive octagonal test section of the Berlin Technical 
College . 

Finally, a new method is presented for testing three- 
dimensional models in wind tunnels with two flexible walls. 

Wall designs are presented by means of which the 
interferences on the model axis can be eliminated. 


*Numbers in margin refer to foreign pagination. 


4 



This work was performed as part of the research project 
"Development of an adaptive wind tunnel test section with 
flexible walls for three-dimensional transsonic flows", 
supported by the Federal Ministry for Research and 
Technology. 


5 


1 


SYMBOLS 


(x >y >Z ) 

( C , n , C ) 

(x,r,e) 

Jvw 

'VW 

( x d’J'd' z d ) 
x >= (x,y,z) 

s .s 1 

B^, Bereich I 
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n, n A 

Iu 

X n 
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Cartesian coordinates 


A 3 

N 

dF 

h 

h(x) 


0 ) 
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B= /1-H* 

U. 


Cylindrical coordinates 
Coordinates of plotted point p 
Coordinates of origin point q 
Coordinates of wall pressure holes d 
Location vector 
Control surface 

Region within the control surfaces 
(near-model region) 

Infinite region outside the control surfaces 
(external field) 

External normal of S 
Internal normal of S 
Beginning of test section 
End of test section 
Initial cross section 

Pinal cross ssction Octagon Test Section 

Mean wall distance from the tunnel 
axis for the j-th surface panel 
Number of panels 
Surface element 
Height of test section 
Wall deflection 

Half-length of the basic wave (design rule) 

Circular frequency 

Incident flow Mach number (uncorrected) 

Compressibility factor 
Incident flow velocity (uncorrected) 

Perturbing velocity potential (normalized 
With y } 
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Perturbing velocity in the x direction 
Perturbing velocity normal to S 
Total axial velocity 
Additional potential (perturbation 
computation) 

Flow function 
Pressure coefficient 
Lift coefficient 
Drtig coefficient 

Pitch moment coefficient at 25% of the 
aerodynamic wing depth 
Ratio of specific heats 
Model angle of incidence 
Compressible source potential 
Source strength 

Influence coefficients (u-component) 
Influence coefficients ( v-component ) 
Influence coefficients for internal flow 
Unit matrix 

Weighting factor in the iterative wall 
adaptation method 
Modified Bessel functions 
Landau symbol 
Tolerance magnitude 


Indices : 

w 


A 

I 

o 


wall-induced 

model-induced, interference-free 
in the direction n A 
in the direction nj 
upper tunnel wall 
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lower tunnel wall 

eff, kor effective (incident flow) 
D designX. 

^ design y' 


design methods 


Strokes above a value indicate averaging or transition to 
conjugated complexes. 
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INTRODUCTION 
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In a wind tunnel, the flow field must satisfy a 
kinematic boundary condition along the test section wall, in 
contrast to a laterally inbounded flow around a model. For 
instance, in test sections with impermeable walls the course 
of the flow lines near the wall is determined solely by the 
wall contour. The flow lines in the wall region, therefore, 
do not generally correspond to those which would be there 
with free flow around the model. This leads to 
interferences in the whole test section flow and introduces 
errors into the aerodynamic data measured with the model. 

The wall-induced perturbation velocity field is normally 
very complex, but it can be at least approximately repres- 
ented as a superimposition of various independent effects. 

The two most important of these are the axial additional 
velocity (blocking) produced by the changed lateral bending 
of the flow lines in the region of the model, and the 
change in model lift due to the diversion of the downflow 
field at the tunnel walls (lift interference) . For example, 
in test sections with flat impermeable walls the blocking 
effect leads to an increase of the effective incident flow 
velocity and, thus, of the stagnation pressure (in the 
subsonic region), while there is generally an increase in 
lift and therefore an increase of the effective model 
incidence angle. The measured data, with interferences, 
then as a first approximation describe the actual 
aerodynamic behavior of the model studied with somewhat 
altered incident flow conditions. Many numerical methods 
have been developed in the course of time for computing the 
corrections and Aa 
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For wind tunnel studies in the transsonic range, where 
the test section walls must be slotted or perforated to 
avoid the transsonic blocking, though, application of 
these correction methods encounters great difficulties. The 
boundary condition for partially permeable tunnel walls can 
only be approximately described mathematically, so that the 
corrections determined on the basis of this information (and 
a representation of the model investigated by means of 
suitable singularities) are likewise troubled with great 
uncertainties . 

Furthermore, one must consider that the "global" 
consideration of the wall interferences in the form of the 
two incident flow parameters and Aa , as described 

so far, is justified only if the wall-induced perturbation 
velocities are nearly constant in the region of the model. 
This is generally the case only with very small blocking 
ratios. On the other hand, in view of ability to transfer 
the model data measured in the wind tunnel to the full-scale 
design, one must certainly try to test the largest possible 
model (Reynolds number similarity) . Because of the 
continuously increasing requirements in recent years for 
quality of wind tunnel measurements in the high subsonic 
region, the design of new test sections which can provide 
model measurements which are free of interferences, or at 
least correctable, with maintenance of a high blocking 
ratio, is becoming more and more pressing. 

This requirement was the cause of the development of the 
so-called "adaptive" wind tunnels, which early in 1970 
realized in design the obvious idea of eliminating 
the cause of the wall interferences by matching the test 
section walls to the laterally unbounded flow around the 
model. First, test sections with two adaptive walls were 
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produced for profile studies, based on highly developed 
computer and control technology. The adaption was produced 
either by deforming the flexible walls into free flight flow 
surfaces or, in the case of ventilated tunnels, by 
appropriate adjustment of the local wall permeability. The 
necessary displacement of the test section walls or the 
adjustment of their porosity must be completely automatic to 
have short test times. For this, it was necessary to solve 
basic design and control technology problems. 

Furthermore, application of this new wind tunnel 
technology would be promising only after a practical method 
had been found for determining the adapted wall shape (see 
Chapter 3). That is, direct calculation of the free flight 
flow surfaces would be very costly because of the normally 
complex model geometry and the friction and compressibility 
effects which occur, and it would not be possible with 
sufficient accuracy, especially in the transsonic speed 
range . 

Testing of the 2D adaptive test sections using some 
standard profiles (NACA 0012, CAST 7 ) categorically 
confirmed that wall adaptation could be done practically, 
thus justifying the added design and electronic cost. On 
the basis of this great success, design of three- 
dimensional adaptive wind tunnels was begun in 1979. Of the 
three designs which have been realized so far, the test 
section with eight flexible walls, constructed at the 
Berlin Technical College, is described in more detail in 
the 4th section. 
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The principal objective of this work, though, is 
presentation of the numerical methods developed at the 
Berlin Technical College for three-dimensional wall adapt- 
ation. Initially, we shall consider generally the 
functional principle of adaptive wind tunnels for 
determining the adapted wall shape. 
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PHYSICAL FUNDAMENTALS OF THE TECHNOLOGY OF ADAPTIVE 
WIND TUNNELS 
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The basic functional principle of adaptive wind tunnels 
was developed by Ferri and Baronti [1] and by Sears [2]. By 
comparing the tunnel flow with the laterally unbounded flow 
around the model they were successful in stating an analytic 
relation by means of which one can determine the wall 
interferences in the test section without explicit knowledge 
of the model being investigated. Derivation of this 
relation and of the iterative wall adaptation process based 
on it is carried out in the following for the general case 
of three-dimensional flows. 

An important point of view here is the division of the 
laterally unbounded flow field (initially seeming rather 
arbitrary) into a region near the model and the so-called 
far- or external; field. In Figure 1 these two flow 
regions, which fuse into each other, are delimited by an 
imaginary control surface. Within the control surface 
(region I) the flow field is very complex for the reasons 
already stated, and can only be approximated with the 
presently known computation methods. In the infinitely 
extended exterior field (region II), though, the model- 
induced perturbations have already decayed sufficiently, at 
sufficiently great lateral distance, that the flow here can 
be considered free of rotation, and the perturbation 
velocity potential can be matched to a good approximation by 
a simple differential equation. Experience has shown that 
the linearized potential equation, or its transsonic 
variant, is sufficient for this. 
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Approximation of the far field by a potential flow now 
leads immediately to an analytical criterion for freedom 
from interference. On the basic of mathematical uniqueness 
rules (see, for instance, [3]) the velocity potential is 
already unambiguously determined if along with the condition 
of vanishing perturbations at infinity, either the model- 
induced perturbation velocity in the incident flow or in the 
normal direction on the control surface S is specified. (S 

can initially be considered as infinitely long.) If one 

\ 

calculates the perturbation velocity potential $ with the 
u-component as a boundary condition on the control surface 
-T7 (S) = u(S), then the matching normal velocity follows 
from differentiation [u] = ypj- (S) and conversely. With 
laterally unbounded flow around the model, the perturbation 
velocity components on the control surface then fulfill a 
compatibility relation which can be evaluated analytically. 

By means of this criterion one can now determine easily 
whether a wind tunnel flow is interference-free, and if 
necessary undertake an adjustment of the adaptive test 
section walls. For this we measure two independent velocity 
components, e. g., the u- and v-flow on the control 
surface which is now assumed to be in the test section. If 
both these measured velocities fulfill the compatibility 
relation above, then the flow around the model in the 
tunnel can be expanded continuously to a laterally unbounded 
flow field, and it is therefore interference-free. 
Conversely, an interference-free tunnel flow can be 
considered part of an unbounded flow field, so that the 
equation derived above applies for the velocity components 
on the control surface. If we designate the measured u- and 
v-flows as Uj(s) and Vj(S), then the flow in the test 
section is exactly interference-free if the relation 
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(3.1) 


V A W] = v i 

applies on the control surface. 

But if the calculated normal velocity v A on the 
control surface, for flow around the model which is 
considered to be extended laterally without limit, is not 
identical with the corresponding measured component Vj, 
v A (S) - Vj(S) ^ 0 , jthen the tunnel flow has interferences 
and correction using the adaptive test section walls is 
necessary. 

The difference, v A (s) - Vj(S) is obviously a measure of 
the magnitude of the wall interferences, and it therefore 
appears reasonable to adjust the contour or the porosity of 
the tunnel walls so as to produce an averaged (with a 
weighting factor K) normal velocity distribution on the 
control surface: 

Vj new = Vj + K(v A - V j) » Ke ( 0 , 1 ) (3.2) 

Now the resulting tangential velocity Uj new (S) is 
measured and another test made by means of an external field 
computation to see whether the flow field within the control 
surface already matches that for laterally unlimited flow 
around the model. The difference, v A [uj new ] - Vj new is 
again used to correct the adaptive test section walls, and 
the entire process is repeated until the wall-induced flows 
have become negligibly small. 

In summary, we have the following iteration scheme for 
wall adaptation. Figure 2: 
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v - v Hi 

V A " v aL u I J 


external fields calculation n=1,2... 


Vi (n*t),»( n ) tK ( VA (n) 


I 


I 


-Vj (n) ), Ke (0,1) (3,3) 


The tunnel flow is considered to be interference-free if the 
change of the v-component in two successive iteration steps, 
and thus the correction of the adjusted wall configuration, 
remains within a specified tolerance e> Q 


(n+1 ) 


- v 


(n) 


<e 


(3.4) 


One should consider that use of the adaptation process 
(3.3) provides no direct information about the model being 
investigated, and is therefore universally useful . The 
complex flow field within the test section is represented 
only by two independent flow quantities, such as the u- and 
v-perturbing velocities, along a control surface near the 
wall. Their repeated measurement along with the so-called 
external field calculation makes possible stepwise reduction 
of the wall interferences. 


It must be emphasized, though, that the convergence 
properties of the control algorithm (3.3) depend strongly on 
the selection of the weighting or control factor K. As a 
result, there is a special practical importance to 
determination of an optimal value K opt which leads to the 
adapted wall configuration with the fewest possible steps. 
Investigations of this sort are performed in Chapter 6 for 
three-dimensional wall adaptations. They show 
simultaneously how (3.3) must be modified with respect to 
rapid convergence (single-step method) . 
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The numerical methods for the external field calculation 
required in each adaptation step are considered extensively 
in the 5th Chapter. In the following section we shall first 
explain the basic problems in the design of three- 
dimensional adaptive wind tunnels. Then the test section 
with eight flexible walls, built at the Berlin Technical 
College, is described in detail. 
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4 . 


3-D ADAPTIVE TEST SECTION AT THE BERLIN TECHNICAL 
COLLEGE 


/13 

In the design of three-dimensional adaptive test 
sections we could utilize the experimental results already 
gained with adaptive profile tunnels. Evaluation of the 
test results was particularly successful with respect to the 
practical feasibility of wall adaptation either with 
flexible impermeable walls or with ventilated test section 
walls . 

The principle of partly permeable walls with locally 
controllable through-flow volumes (controlled either by 
variation of the wall porosity or by individual matching of 
the counterpressure in the sectionally divided plenum 
chambers) has been applied particularly in the USA. It 
appeared advantageous, especially economically, because the 
usual transsonic test sections are slotted or perforated, 
and could therefore be re-equipped relatively easily. 

In comparison to the adaptive profile tunnels with 
flexible walls, built in Europe, this technique showed 
significant practical disadvantages from the beginning. 

This became particularly apparent in the measurement of the 
near-wall perturbation velocity distribution at both sides of 
the profile in the flow, required at each adaptation step. 

With impermeable walls, the two desired velocity components 
can be determined directly from the wall shape (v = h') and 
the pressure distribution prevail ing there (u = - Cp/2). 

But with ventilated tunnels, measurements within the flow 
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field are necessary because of the severly inhomogeneous 
flow near the walls. These can, to be sure, be performed 
without contact by means of lasers, but this method is 
relatively time-consuming and is defensible only for 
adaptive profile tunnels. 

The problem of rapid and precise measurement becomes 
more important in design of a three-dimensional adaptive 
test section because of the increased amount of data 
(longitudinal and normal perturbing velocity distribution 
along a control surface surrounding the model). It was 
principally for this reason that in Germany (at the Berlin 
Technical College and later at DFVLR, Gflttingen) the 
decision was for design on the basis of impermeable, 
deformable walls [4]. (Other reasons were the higher power 
consumption and the higher noise level of ventilated test 
sections. ) 

At about the same time in the USA, though, work was 
started on construction of a three-dimensional adaptive 
tunnel with four ventilated walls [5]. A new kind of 
measuring system, which has become known as the "Calspan 
Pipe" was developed to determine the perturbing velocity 
components. It consists essentially of two metal tubes with 
static pressure holes, which can be turned about the 
longitudinal axis of the wind tunnel model. Figure 3. To be 
sure, experience so far has not been satisfactory. That, 
along with other problems, has as yet prevented a 
convincing functional demonstration of this test section 
C6]. 


In contrast, it has been possible to carry out initial 
successful wall adaptations with the two other 3D adaptive 
test sections, and to demonstrate the general applicability 
of the technique, which has in the meantime been established 
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for profile studies, for three-dimensional flow cases as 
well (see Chapter 9). The so-called "extensible adaptive 
test sections" built by the DFVLR in Gottingen is 
described in [7] along with the most important experimental 
results. In the following we shall consider the design 
details of the test system developed at the Berlin Technical 
College. 

Design of the test section began with the basic idea 
that an ideal three-dimensional wall adaptation is hardly 
achievable in design, and not achievable at all at least 
for subsonic flows. Important viewpoints for the design, 
then, were the provision of the least possible mechanical 
cost while simultaneously providing adequate deformability 
of the test section walls. This compromise which was 
attempted led finally to design of the "octagon" test 
section with eight flexible walls [8]. 

Figure 4 shows the octagonal test section cross section 
formed from the eight walls and the layers between them 
(the model is described in Chapter 9). The principal 
dimensions (height 15 cm, width 18 cm) were chosen in the 
ratio of Is 1.2 with respect to a planned full-scale 
"European Transsonic Wind Tunnel", ETW. The test section 
length of 83 cm arises from the dimensions of the 
connections to the Technical College wind tunnel . Each 
flexible wall is individually deformable by means of 10 
positioning elements in each (9 in the upper and lower 
walls), positioned longitudinally. The deformation and 
position control are fully automatic C8]. Figure 5 is an 
outer view of tine test section showing the errengeinent of 
the positioning motors. The wall pressure distribution is 
measured by a total of 192 (6x 25; 2 x 21) pressure holes 
along the center lines of the flexible walls. 
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It was practical to choose the flow tube formed by the 
straight, undeformed test section walls as the control 
surface for the external field calculation (see Chapter 3). 
(Selection of a fixed interation-independent control surface 
is important for performing the adaptation quickly.) The 
perturbation velocity components needed there can be taken 
unchanged from the wall measurements because of the expected 
small contour changes. (The effect of the wall boundary 
layer is considered approximately by calibration of the 
empty test section. Chapter 9.2.) 

The next sections describe extensively the numerical 
methods developed for the external field calculation. 
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5. 


NUMERICAL PEFORMANCE OF THE EXTERNAL FIELD 
COMPUTATION 


5.1 Basic Comments 


The boundary value problem to be solved in the external 
field computation is, using the "theory of small 
disturbances" (Chapter 3), Figure 6: 


+ *yy * = - 


<t> (X > ) * / u T (t ,y ,z )dt 


0 

(k .1 ) M“* x * xx 


* 


— > 
X 


eS 


subsonic 

transsonic 


(a) 

(b) 


(5.1) 


♦ = o (- 4 -), 
1 5* I 



- > 


oo 
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We seek the normal perturbation velocity components 

v a[ u j] = on tbe control surface S. Clearly explained, 
this is a matter of calculating that wall shape which would 
experience the measured pressure distribution in the fictive 
external flow. 

Then the (unfavorable) physical properties of this 
three-dimensional design problem are demonstrated by means 
of a singularity representation for the solution in subsonic 
flow. In the next chapter we then present the design method 
supported by computation developed for this purpose. 
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We generate the desired potential by a source-sink 
distribution over the control surface. (Sources and sinks 
are preferred to a dipole coverage because of the lesser 
"leakage” [343). Then we get the unambiguously solvable 
integral equation of the first type [33/ 

i 

- *7$ 1 a (q)g(p>q) dc i = *(p)> p^s ( 5 . 2 ) 

s 

for which the kernel g(p#q) of the (compressible) source 
potential corresponds to 

g ( p *q ) = { p- < x p - x q ) 2 + (y p 

+ (z p " ^ = 1 

Using the relation (5.2), we can demonstrate the conjecture 
already stated in [103, that the three-dimensional design 
problem is ill-posed for this particular case, by means of 
a simple perturbation computation: 

In order to be able to investigate conveniently the 
effect of small pressure fluctuations on the resulting wall 
shape, we approximate the control surface. Figure 6, by a 
circular cylinder, and limit ourselves to local 
rotationally symmetric perturbations in the basic Uj 
component. The induced supplemental potential ^ is then 
constant downstream from the perturbed site. In this 
region, according to Equation (5.2) it produces a likewise 
(nearly) constant source distribution. (Differentiation of 
(5.2) and then partial integration leads to the relation 

" TTp £ -^7^g(p»q)dq = Uj(p) , peS (5.4) 

b q 


y ) 2 
y q 


M 2 

oo 


(5.3) 
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and accordingly ° = constant for Uj = 0.) its 

induction effect in the external field can be expressed by a 
line source on the center axis of the control surface 
running from x * -» to x = +« . Then, from * Tif ? nr 

which is linear at the location of the control surface 
one obtains for the perturbation-induced v-component 

v A («+) = «h'(x) a = const,' <5.5) 

As conjectured, a local pressure perturbation finally leads 
to large global contour changes, (Th • 


This situation is particularly important for practically 
performing the external field computation because the 
control surface potential $(p) , because the Uj 

distribution is known only at the pressure holes, can be 
determined only up to one constant: 


<f»(p) 



1 u i ^ t »y p » z p^ dt+ 


1 u i (t -v z P )dt ’ 

X, 


peS 


(5.6) 


x A : beginning of test section 

The initial potential <f> (x^»Yp which must be 

determined by extrapolation can exert a considerable effect 
on the desired wall shape. With the C5 model, for example, 
there can be differences of 30% even in the interference- 
free case in the wall deflections at the end of the 
test section, depending on whether or not its flow effect 
is considered (<f>(x^,R) = 0) . (See Chapter 8.2.) 

Furthermore, the numerical integration of the Uj 
distribution on the basis of the relatively small number of 
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support points (25 or 21 pressure holes per flexible wall, 
see Section 4) leads to an error (which grows in the 
direction of flow) in the potential determination, and 
therefore leads to an incorrect wall contour, especially at 
the end of the test section. 

Direct use of the (extrapolated) Uj component as a 
boundary condition, on the other hand, leads to an 
integral equation which cannot be solved unambiguously: 

1 TxT 9(P»9) dc l = Uj(p), p e s (5.7) 

S P 

i 

and therefore, after discretizing, to a singular equation 
system for the source intensity ° (see (5.4)). 

Nevertheless, in order to be able to test the computed 
wall deflections in the simplest possible way in the present 
application case, the design problem is once more solved 
iteratively using a design method . Here the desired wall 
shape is determined by stepwise modification of a starting 
contour. The necessary geometry changes are determined from 
the deviation of the actual pressure distribution induced by 
the particular existing contour in the external field from 
the specified final values Uj. The design rule developed 
for this is derived in the next chapter. 

To compute the wall pressure distributions one can 
select either the singularities method or the TSP method. 
This computation problem in the subsonic range leads to a 
numerically advantageous integral equation of the second 
type (Section 5.2.2). Even in the transsonic range it is 
significantly easier to solve than the comparable design 
problem [5,11] (Chapter 5.2.3.) 
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In the subsonic range the initial contour can be the 
straight, undeflected test section wall. In using the TSP 
method, though, it is better with respect to computing time 
to provide a good starting value using the panel program. 

/21 

5 . 2 Design Process 


5.2.1 Derivation of the Design Rule 

In Figure 7 we show the outer field calculation v^u^] carried 
out using the design process. In this chapter we will give a 
detailed description of the design rule. 


The basic idea in the derivation was the observation that the 
pressure distribution along the (slightly) deflected test 
section walls corresponds approximately to that along the outer 
side of a geometrically similarly deformed circular cylinder. 

If we approximate the wall contour (averaged in the circumferential 
direction) by means of a sine half wave 

MX) = csin (£x) , J „ sL , c = const ’ (5 ' 8) 

then for subsonic flow (^<1) we obtain a simple analytical relation- 
ship between the flow angle change VjJ = h" and the axial perturbation 
speed u D , Figure 8: 


First, for the potential distribution in the outer field we 
have [12] 


4> D (x,r) 


C K 0 ( r 8r > 

8 7^7 cos(fx) ' r * * 


(5.9) 


and from this after differentiation at the position r = R and 
then solving for v^ we find 
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eirK 'i ( r bR) 

L K q (£ PR) 


u D (x , R) 


(5.10) 


(5.10) corresponds already essentially to the desired design rule 
for the wall contour calculation, but the wavelength L has to be 
matched to the prevailing test conditions. 

In general the measured wall pressure distribution (and also 
the desired wall contour) cannot be exactly defined by a single 
sine halfwave. Instead it will contain high frequency contri- 
butions vto n = n L"* n ~ 2 S 3,... in the sense of a Fourier series 
expansion. In order to nevertheless be able to obtain a proportionality 
between v^ and u D , the influence of L on the proportionality factor 
in (5.10) was investigated in detail. The results are given for 

P 0.7 as an example and a "basic wavelength" L = 1000 in 

the following table which is matched to the length of the octagon 
test section (half wavelength) „ 


L 

z = { eR 

K 1 (Z) 

e, M z) 

yw 

-ri^TiT 

125 

1.492 

1.29 

- 0.023 

250 

0.746 

1.55 

- 0.014 

500 

0.373 

2.04 

- 0.009 

1000 

0.186 

.L 

2.89 

- 0.007 

— — 


As can be seen, the proportionality factor increases in /23 

magnitude for the harmonics but the overall change remains 
relatively small. Considering the dominance of the basic wave 
at least for pure displacement flows and in order to provide 
monotonic convergence of the design cycle, we decided to use the 
design rule given above with a fixed value L = 1000 (see Chapter 8.2 
for a discussion of convergence behavior) . 
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Because of 


irgR 


< 0.3 for L = 1000 the Bessel functions in 


(5.10) can be expressed by their symptotic expansions for a small 
argument [13]. The desired design rule therefore has the following 
final form 


vj(x) = 


1 


^ ln ^ 1 000 ^ “ M ~ 21 ^ 


u d (x,R), M^<1 


(5.11) 


This corresponds essentially to the differentiated form of (5.5) 
and therefore not only validates this approximation formula but also 
the discussion there about the "ill posedness" of the design problem 
(5.1) . The constant which occurs in the integration of 


v D (x) = / Vq( x) dx + C 1 


(5.12) 


can be ignored for sufficient upstream extrapolation. 

5.2.2 Post-computation process for the subsonic region (Panel method) 

The calculation of the axial perturbation speed along the wall 
contours obtained during the individual design steps can be reduced 
to the following Neumann boundary value problem. /24 

02 *xx + fyy + *zz = 0 » e 2 = 1 - M* >0 

f-fj- (x > ) = fiv D (x > ), XeS 

j (5.13) 

= 0 ^ | > | ^ , | x > | — > «° 

where n j is the outer normal along the control surface S (Figure 6) . 

If we produce the desired potential <{> by a source-sink distribution 
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which is continuous on S (Formula 5.2, 5.3), then we obtain the 
source density o as the solution of the following integral equation 
of the second kind 


1 °(p) - f °(q)|n g(p><?)dg: = «v D (p), peS 

S lb . 14 

. o 

The quantity y describes the induction effect of the source 
intensity o(p)dF located at the target point which is not con- 
tained in the integral. 


As an illustration, consider a source occupation in the plane z^ =0 
with a target point contained in it ( x p»Vp»0) • Then we have 

9(p,q) = (x p - V 2 + (y p - y q ) 2 >' * 

and therefore 


/ o(q)|=“ g(p,q)dq = °- 

3Z p 

Z q =0 


The integral, even though it expands over the entire plane, 
therefore only gives the influence (zero in this case) of the source 
distribution outside of the target point. The local part has to 
be determined separately. Application of the Gauss theorem [3], /25 

on the charge element located in p, then immediately gives the 
relationship 


/ 3 <t> 
v an, 


(P) 


InY <»> 


dF 


a(p)dF 


and therefore 



- y 0 ( P ) ) 


With this extra term, the linear equation system which approxi- 
mates (5.14) takes on a very favorable numerical property (see 
below) . 
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In order to perform the discretization, the control surface 
is divided up into flat quadratic surface elements , the so-called 
panels 


3 ■ U 

j = 1 


4 o 


(5.15) 


and the desired singularity intensity is approximated by a step 
function (panel method of first order) 


(q) 


o . s const. , 
J 


qeAj , j = 1 N 


(5.16) 


For the normal speed V induced in PeS we then find 

1 ci 


V ( P ; Oy, 


'N 


) = Y o(p) 


1 

4 it B 


N 


0 = 1 


. ; g(p.q) d( i 


(5.17) 


A o 


The constants 
boundary conditions 


have to be matched so that the 
are approximated as well as possible. 
This is done most simply by the so-called "point-matching", [3], 
in which the approximation function is determined in such a way 
that it exactly agrees with the prescribed distribution at fixed 
points of the definition region. In the present case one selects 
a "control point" p^ on each panel and then from the requirement 
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v ( P -j > o ^ » • • • »Oj|) - ^Vp(P-j) > 1 — 1» N 


(5.18) 


one finds the linear equation system 

j, A ij ’j ■ 6 , D.i . 


with the influence coefficient 


i J 



i = 1 , ...» N 


1 = J 


4tt B 


/ g(p i ,q)dq i 4 J 


A j 


(5.19) 


(5.20) 
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The calculation of A^ ^ is performed in the Appendix (Section 13) . 

In addition to the exact formulas, we give simple approximations 
in the form of a multipole expansion, which can be used for a large 
distance between the target point and the induced panel, [17]. 

. 1 /2 

The diagonal elements ■jr are much larger compared 

with the other coefficients and are very advantageous for the 
numerical solution of (5.19). Therefore iteration methods can be 
used for this. This is important for practical applications, because 
due to the memory requirements of direct methods (Gauss-Elimination, 
etc), there can be restrictions on the number of surface segments. 

(The coefficient matrices which occur during the discretisation of 

the integral equations of the first kind, see for example (5.2) , 

have as a rule no dominating diagonal and can only be solved directly) . 

However, we should not forget the fact, that the values o^, ...» 
calculated from (5.19) only represent a usable approximation for the 
desired force intensity o (Equation (5.14) if the panel coverage 
of the control surface is appropriate). This discretisation error 


<5 Q = / | o ( q ) - o(q) | dq 

S 


(5.21) 


can be simply reduced by increasing the number of elements. But 
especially in the present application, restrictions can be 
brought about by the limited performance of the wind tunnel computer. 
However, the deflections of the test section walls only have to be 
determined within the adjustment accuracy of the past transducers 
( maximum | ) (it is 0.03 mm in the model area and 0.07 mm for the 
octagon test section) . (The influence of the wall contour accuracies 
on the model flow is discussed in Section nine) . 


In order to find a panel coverage which is efficient for this, 
we first determine the adapter wall shape for two typical flow 
cases (C5 and F4 model represented by singularities) . 
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This was done with a very fine segmentation of the control surface, /28 
(16 x 30 elements) using a large computer. After this the number of 
panels was reduced with consideration of the local source gradient 
according to the tolerance specified above. The optimized paneling 
is given in Figure 9. It consists of 16 strips (according to the 
8 flexible walls with intermediate lamella pieces (Figure 6) , which 
are 15 elements each of different lengths. The detailed discretisation 
in the circumferential direction is required because of the 
necessity for covering the vortex train which induces lift along the 
rear test section region. However, since only symmetric models with- 
out slip angles are investigated, the resulting linear equation 
system (5.19 can be reduced to only 9 x 15 - 135 nonredundant individual 
equations . 


The numerical solution of the boundary value problem (5.13) 
is 

4>(p) = ~ > peSuB e (5.22) 

S 

where * is the calculated approximation according to (5.16) and 
(5.19). The desired perturbation speed 6d e along the control 
surface is then obtained by differentiation 

6 u E (p) = ^ = " 47 b 9(p»q)dq, peS (5.23) 

S P 

An evaluation along the control points P^* 1 = 1 > •••» ” 
gives 


6 u 


E,i 


N 

z 

j=1 



(5.24) 


with the influence coefficient 



1 

1 


3 



9 ( P j » Q ) d q 
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(5.25) 
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and the calculation is also given in the Appendix. 


A combination of (5.19) and (5.24) then leads to the 
very simple linear relationship 

6u^ = BA _1 6Vp 


(5.26) 


where the matrix BA - "*" formed from the influence coefficients 
can be determined before the test. The post calculation task in 
the subsonic case therefore can be reduced to the use of a panel 
method and a simple matrix vector multiplication. 


5.2.3 Post-computation process for the subsonic region (TSP method) 


•In order to also discuss the case of incident flow near the 
speed of sound, the post calculation part of the design cycle 
can also be done using a TSP method. The corresponding computer pro- 
gram was made available by AEDC | 5 | and will be briefly discussed. 


We solved th6 classical Guderley-von Karman equation [10] in the 
conservative formulation 


[<i-»L>*x 



with the boundary condition 

(x, R, 0) = sv D (x,e) 


/30 


and the circular cylindrical assumed control surface 

(Radius R) and 

4 ->0 r -> 00 


-> 0 1 x 1 -> - 


(5.27) 


in the far field. 
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Based on the symmetry relationships ($(x,r,e) = <|>(x,r,-0)) 
we finally find 

<|) 0 (x,r,O) = 4> e ( x ,r ,it) = 0 , r > R 

The difference network was matched to the present application 
case and now extends in the range 

1 0<x <940 1 x 1 83<r* 332 1 x |0 <O<tt| with 31 x 10 x 11 = 3410 

points. The discretisation of the space derivatives and the 
solution of the resulting difference equations is done using standard 
technique (finite difference method of Murman, SLOR- relaxation 
method) which are discussed in detail in [10] and [18], 

The numerical complexity for calculating the desired perturbation 
speed distribution 6u^(S) = (S) is of course much more than 

for the subsonic case. Therefore, considering the desired Online 
wall adaptation, one should use the simple relationship (5.26) as 
much as possible. 

The range of the validity of the linear potential equation is dis- 
cussed in detail in Section 9. 

6. OPTIMIZATION OF THE RATE OF CONVERGENCE FOR THE WALL /31 

ADAPTATION PROCESS 


6.1 Subsonic Case 


6.1.1. Convergence Behavior for Constant Regulation Factor 

By a suitable selection of the regulation factor K we wish to 
increase the convergence rate of iteration wall adaptation method 
(3.3) similar to the two-dimensional case. 
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A first indication about the effect of K on the adaptation 
method was found from demonstration tests in the octagon test 
section of the TU Berlin. It was found that an adaptation of the 
test section walls was only possible with very small regulation 
factors (K < 0.10) and relatively numerous (*10) iteration 
steps. For larger values of K one found oscillating wall contours 
with rapidly increasing deflections, which already lie outside of 
the permissible displacement paths of the past transducers after a 
few iteration steps. 


In order to make a systematic investigation of this experimental 
finding, we performed a numerical simulation of the regulation 
algorithm (Figure 2). The wall pressure distribution u^. ^ n+1 ^ i s 
not determined by means of the wind tunnel tests in the individual 
adaptation steps, but instead using a natural assumption about the 
given wall contour v.^ . 


Let us start the simulation with the measured values 
and define 


u j ^ = u i ^ ^ + j ^ 


6V 


I 


( 2 ) 




(6.1) /32 


( 2 ) ( 2 ) 

Then one can establish between 6u^ and Sv^ (accordingly for 

6Uj n+ ^ and Sv^. in the following iterations) an analytical 

relationship corresponding to (5.26), if we assume that the flow 

( 2 ) 

field induced by contour change Sv-j. is irrotational and the 
potential ^(2) satisfies differential equation (5.1a). These 
conditions are satisfied because of the assumed small wall deflections, 
to a great degree of accuracy [7,9]. 

The trial solution (the loop index is omitted) 


6<<,( p) = - ^ / C(q)g(p,q)dq , p eB .US 
S 1 


( 6 . 2 ) 
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then leads to the following integral equation of the second kind for 

C 

1 ‘ 77? 1 ?jr~ g(p»q)dq = - 6vj(p), peS (6.3) 

s 1 

when we finally find the desired increment Su^. as 

6U I (P) = " 77$ ; 3T" g(p,q)dq, peS (6.4) 

S P 


After performing the discretisation discussed in Chapter (5.2.2) 
we then find the two linear relationships 


A , C > 


6uJ 


= ' 6V I 


= Bc > 


(6.5) 


and the influence coefficients Ft . can be expressed by the 


already defined (5.20) 


A 'u 


L 


A u 


- A U 


i = J 
i + J 


( 6 . 6 ) 


A comparison of (6.5) with the corresponding relationships 
(5.19) and (5.24) for the outer flow 
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Ao” 


6U A 


**A 


Ba* 


(6.7) 


( 6 V 




+ ) 


immediately shows that the same wall contour change VUV A 
in the outer field produces a much lower pressure fluctuation 
than inside the test section. This then makes understandable 
the adaptation behavior discussed above. 


This state of affairs is found from the estimation 


A i j > 0 * i.J = 1 , . . . , N 


( 6 . 8 ) 
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which follows immediately from the convex nature of the volume 
inside the control surface. Figure 9. In the outer field 
calculation (6.7) the singularities therefore mutually support one 
another with their induction effect. This means that for producing 
prescribed deflection 6v^ , a smaller charge density a is 

required. However, for the inner flow we find 


N 


I 

j=1 


1 J 


i 


N 


(6.9) 


a 


Equation (6.3) for $ v j(p) = 0 has a nontrivial solution ^ 

£=const because a constant source distribution produces a constant 
potential along the control surface and therefore inside of it as well. 
See (5.4)). In this case the induced normal speeds approximately 
cancel and one requires a large source-sink intensity in order 
to obtain the same contour change +6v^. Since 5u^ and 5u^ are deter- 
mined from the same matrix B, we obtain the assumed physical state 
of affairs. 


The results of the simulation are given in Figure 10 for a 
typical flow case [19], (In the Appendix we also give an analytical 
investigation of the convergent behavior) . As expected, the adaptation 
method only converges for very small regulation factors K ^ 0.06 
One remarkable feature is the large increase in the number of 
iterations for factors larger than the optimum value ^opt * 
and the ensuing transition to divergence. Report [9] gives an 
account of a similar result during the adaptation of three-dimensional 
flows around wings. The numerical simulation therefore not only 
confirms the experimental results but also shows that even if one 
uses an optimum regulation factor, the matching of the test section 
walls requires an excessive number of intermediate steps and there- 
fore long wind tunnel testing times. 


In the next section we will show that by means of a differentiated 
weighting of the channel line flow and the outer flow, one can 
achieve a drastic increase in the convergence rate. 
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6 . 1.2 


Single-Step Method 
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We again start with the measured values Uj (S) , (S) and the 

outer field calculation • We then investigate how one can 

obtain the best possible approximation for the desired adapted 
wall shape. 

For the mathematical description of the problem we will 
formally decompose u^ and v^ into the model part without inter- 
ference and the influence of the wall. 


u T = u_ + u., 


v t = v + v.. 
l «> w 


( 6 . 10 ) 


Based on the linearity of the outer flow we then find 


v A[ u l] V I = v a> + V A^ U W-1 " v » “ V W = v a£ U kJ ' V W (6.11) 

The design task v a[ u W-J can then be solved in two steps 
according to (6.7). In this case the singularity “ intensity o 
is of course unknown. On the other hand, the wall-induced pertur- 
bation speed field satisfies the linear potential equation (5.1. a), 
to a good approximation, so that we can also use the relationships 
(6.5) discussed in Section 6.1.1. 


Overall, we obtain the following four equations 


Ao = v aK|] A'c = - v w 

Bo ’ u i; Be =, u w 


( 6 . 12 ) 
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and from this because of (5.4), (6.9), (6.6), (5.20) 
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C = a + c , c = const. 


(6.13) 


Va[u w ] - V W = A a + A'(a+c) * (A + A')o = Ea = o (6.14) 
where E is the unit matrix. 

Together with (6.11) we then obtain the following formula 
for the wall dependent perturbation potential 

♦u(P> ‘ - TIs ‘ tV A (q) ’ V I (q)) 9(P ’ q) dq (6 . 15) 

3 

(for V A^ U I^ we will write in the following. as an 

abbreviation) . 

The (known) difference V A^) “ v j(^) therefore is not 
only a "measure" for the wall interference (see (3.2)), but 
also determines it uniquely. It should be noted that this statement 
does not require any additional assumptions about the model- 
induced perturbation field. 

Application of (6.14) directly leads to the 

Single step form 

= v j - Vy. = Vj + A' ( V^-Vj ) = Vj + (E-A) j) (6.16) 

The adapted wall contour can therefore be determined already 
from the measured values of a single wind tunnel test (within the 
range of validity of the linear potential equation) . /37 

This can be done even without numerical additional complexity 
because the matrix A already has to be calculated once (Chapter 
5.2.2). As a comparison which (3.3) shows, the enormous gain in 
convergence is based on the replacement of the constant regulation 
factor by the "regulation matrix" E-A which considers the inci- 
dent Mach number and the control surface geometry. 
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In this connection we should emphasize that already before a 
single step method for three-dimensional wall adaptation was 
developed at the DFVLR and was used successfully [7], The 
Fourier trial solution selected there for representing the wall- 
induced perturbation potential however, has substantial practical 
disadvantages compared with the singularity method, such as for 
example a complicated formalism and the fact that one cannot 
transfer to other control surface cross sections. A comparison of 
the two methods using measured values from the "expandable 
adaptive test sections" of the DFVLR is performed in Section 8.3. 
After this, we report about experimental results in the octagon 
test section. 

Since relationship (6.15) can be used without model representa- 
tion, it is very useful even alone. One can therefore develop a 
high performance correction method. In addition, it is the 
foundation of the new method for testing 3D bodies in wind tunnels 
with two flexible walls, discussed in Chapter 10. 

6.2 Transsonic Case 


The regulation matrix E - A defined in the previous section 
depends continuously on the intimate Mach number and for — > 1 

goes to a well defined limit value (formerly discussed in Appendix) . 
Since the TSP method (Chapter 5.2.3) only can be used for /38 

flow near the speed of sound, considering (6.16) the following 
algorithm of wall adaptation promises a fast conversion. 


Vj (n+1) = Vj (n) + ( E-A) (v A [Uj (n) ] - Vj (n) ), 

M „ =1 (6.17) 


n = 1 , 2 , 


In contrast to the subsonic case, more iteration steps will be 
necessary in general. 
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7. EXTRAPOLATION AND CORRECTION OF THE MEASURED PRESSURE DISTRIBUTION 


In order to transfer the Cp distribution determined along the 
wall pressure traps to the panel control points, Figure 9, 
the measured values have to be interpolated in two dimensions 
and have to be extrapolated beyond the test section region.. By 
a continuation upstream at the same time, we wish to specify the 
effective incident speed U e f£ and if necessary a correction to the 
boundary conditions u^ have to be performed. 

In the derivation of the extrapolation method developed for 
the octagon test section, we started with the usual singularity 
representation of the model and the flow [21] : 

Point Dipole for Simulation of Model Displacement 

Horseshoe Vortex for Representing the Model Lift 

Point Source for Representing the Wake (Including 

Model Spear) 

Because of the different induction behavior of these singulari- 
ties, it is immediately found that the perturbation pressure distri- 
bution which prevailed at both test section ends can be attributed 
mainly to the influence of the wake. In the extrapolation of the 
wall measured values, we only have to consider the source effect. 
(The vortices in the wake are assumed to align parallel to the 
incident flow. Transverse flow components are ignored) . This pro- 
cedure not only gives physically reasonable results, but also 
avoids the numerical difficulties associated with a detailed flow 
simulation [22] . 

The calculation of the wake source has no problems because due 
to the impermeable channel walls, it is given independent of their 
deflection, by the fictitious mass flux difference between the 
test section beginning (x=x^) and the final cross section (x=x £ ) . 
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If we call the axial total speed distribution measured at the 
pressure caps d^ j , i = 1,...,8, , j=l,...,26 U g es ^ij^» then 

by introducing the nominal (uncorrected) incident speeds (Mach number 
MJ we find 


U 


U ges^ d i 1 ^ 


(7.1) 


for the (model and wall-induced) perturbation part 


U I (7.2) 

i 

By using the Gauss theorem we find [3] the following for the 
desired force intensity 


Uj ' (x A ,y,z)df 
(x D,26^ * F E 

(7.3) 

and this can be calculated without any problems during test operations . 

The induction effect of this source of course depends on the 
contour of the test section walls and therefore we will separately 
investigate the two cases of strong and weak wall influence : 

(aerodynamically) straight walls: /41 

The source produces a constant perturbation speed which is the 
same (in magnitude) over the cross section at the two test section 
ends due to symmetry, 6u. Together with (7.3) we find 

6u(x E ) = - 6u(x a ) = -iy Uj'tXp 26 ) (7.4) 


f Uj (x E ,y ,z ) df - j 


_ 1 
~ £ 


8 

i = 1 


U I ' (d i , 26 ) 
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In this wall configuration, consequently we find a non- 
vanishing full effect of the model in the flow, which has to be 
considered by a corresponding correction of the incident flow 
conditions. Because of (7.1) and (7.4) the corrected (effective) 
incident speed U is then calculated as 

U eff ' U » - 5 u < x A> * U. ♦ 7 S l‘(K D>26) : (7.5) 

which brings about a change in the incident flow Mach number of 

. (1 , M„ 2 ) H. • \ V (x D.26> (7.6) 


[ 21 ] . 


From (7.5) and (7.2) we finally obtain the corrected axial 
perturbation speed 

U I (d U ) = U ges (d 1j ) " U eff = U I * < d i j > - 7 5 I ' (x D.26 ) (7 * 7) 

which according to (7.4) can be continued beyond the test section 
region using 


Uj(x) = i 


7 U I ' ^ X D , 26 ^ * 


7 U I ' (x D,2G^ * 


x < x 


x > x 


D,1 


D , 26 


(7.8) 
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The transfer of this extrapolated Uj distribution to the panel 
control points is then done using a two-dimensional cubic spline 
[23] . 


Deflected walls: 


Based on the single step formerly discussed in Section 6.1.2, 
we will restrict ourselves to the case of almost adapted test 
section walls. 
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We obtain the following for the perturbation speed induced 
upstream by the weight source 

6 u(x a)~ TvV 0.9) 

where (x = 600, 0, 0) is the position of the singularity in the 
NL 

tunnel . 

A connection of the incident flow conditions therefore is in 
general not necessary if there is only a weak wall influence. Almost 
the entire source induced mass flux (7.3) therefore passes the 
rear test section cross section, so that the following extrapolation 
of the perturbation speed distribution (7.2) makes physical sense. 


Uj ' (x) 


i 


‘D,1 


(7.10) 


x " X D , 26 


The interpolation of the boundary conditions to the panel control 
points is again done using the already mentioned 2D-spline method. 
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8. TESTING OF THE ADAPTATION SOFTWARE 


8.1 Structure of the Program Systems 

Figure 11 shows the programming system for the three-dimensional 
wall adaptation in the subsonic case. It consists of three separate 
main programs, which are coupled with data transfer. 

The geometry program "GEO" defines the panelling of the 
control surface using the reading surface points (Figure 9). The 
calculated coordinates of the panel corners are transferred together 
with the selected incident Mach number using the program "EINFL" 
for determining the influence coefficients A.., B. . (Chapter 5.2.2). 

The storage of the calculated values is done in the form of the two 
matrices BA ^ and A considering their further use in the Program 
"ADAPl" . The structure of ADAPl is given in Figure 12. The three 
important features are: 

(1) Extrapolation and correction of the measured pressure 
distribution (Chapter 7) 

Transfer of u^ and v^. to the panel control points 

(2) Outerfield calculation v A [u T l using the design method 
(Chapter 5.2). 

(3) The evaluation of the single step formula (Chapter 6.1.2). 

With the calculation of the two influence matrices the program 
Parts (2) and (3) are reduced to a simple matrix vector multiplication, 
so that a fast calculation of the adapted wall contour is assured. 

(see Chapter 9) . 

For the adaptation in the transonnic case, a separate computer 
program "ADAP2" was written. Its structure is similar to that 
of "ADAPl", but within the framework of the outerfield calculation 
we used the difference method developed by AEDC for the transsonic 
case (Chapter5 . 2 . 3) . The regulation matrix (E - A) ^ =1 (Chapter 

oo 

6.2) is determined in another sub-program (as can be seen from the 
formulas in the Appendix, this can be reduced to the size 5 x 5 for 
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However, in this case, the wall contour calculation requires 
a large amount of calculation time and memory and therefore cannot 
be done on line (Chapter 9). 

8 . 2 Analytical Test Case 

Ths computer programs for wall adaptation are tested using 
analytical test cases and by close calculation of known experimental 
results. The calculations used for control of the TSP method are 
discussed in detail in [29]. In this section we will discuss the 
theoretical test nans with "ADAPl" . 


The simulation of an interference flow field in the octagon 
test section was performed using singularities. As a model we used 
the rotationally symmetric ONERA C5 body. Figure 16a, which has 
displacement effects which can be represented in a known way by 
sources and sinks along its axis of rotation [14]. We then 
find the following closed solution for the model induced potential. 


<J>(x,r) = 


77 


FLU) 
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(x-£) 2 + s 2 r 2 I * 


d5 


B 2 = 1 -M 


( 8 . 1 ) 


where F ' (£) is the cross section area change and l is the length 

of the body. If we introduce 


F'U) = \ irdUJd'U) (8.21 


the model thickness distribution d(£) as a geometric quantity in 
(8.1), then we find from it after partial integration the easy to 
evaluate formula 


<j>(x,r) 


1 

TE 


a 

I 

o 


d 2 U) (x-c)^ 

[(x-s) 2 4B 2 r 2 ] T 


d5 


1 d 2 U) 

^ [( x -£) 2 + B 2 r 2 ]^ 


(8.3) 
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The term outside of the integral occurs because the 
rear model diameter (d(£) was not made zero in order to simulate 
that model sphere. 

From the two velocity components obtained by differentiation / 4 ^ 

of (8.3) we can calculate as usual the static pressure and the 
interference -free streamline course everywhere in the flow. 

In order to represent the influence of the wall, we used the 
panel method discussed in Chapter 5.2.2. The singularity intensity 
along the panels was determined in such a way using Equation (6.3) so 
that after superposition with the model induced flow field (8.3), 
the kinematic boundary condition 


v 




(8.4) 


is exactly satisfied along the control surface (C5-body in the 
octagon test section with non-def lected walls) . 


The resulting u-distribution at the points of the wall pressure 
shafts 


ii _ 9<J> . 3 tf’W 
i ax ax 


(8.5) 


together with the incident flow Mach number M and the values 

00 

v-j. = 0 (Equation 8.4) were then the input data for the adaptation 
program "ADAPl" . The adapted wall shape calculated with (single 
step) method in addition to the corresponding pressure distribution 
was then compared with the exact interference free variables calcu- 
lated from (8.3) as a check. 


Figure 13 gives results for a test section wall at = 0.7. 

The agreement of the calculated wall pressure distributions is 
satisfactory overall, even though the extreme values were not 
exactly achieved. The variation of the wall deflections calculated 
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using " ADAP 1 " agree also quite well with the corresponding adapted 
contour. The somewhat larger deviations at the end of the 

test section are still acceptable considering the past transducer 
tolerance of 0.07 mm there. Also, one can achieve an improved /48 

reproduction by an exact determination of the Uj-distribution 
input to "ADAPl" (8.5) (more panels for representing the influence 
of the wall) . 

Figure 14 shows the typical convergence behavior of the 
design method represented using values for the analytical flow 
case discussed above (M^ = 0.7). The accuracy e is defined by the 
mean relative deviation between the actual and nominal distribu- 
tion at the individual designed steps (see Figure 7) 

, N |U,(k)-U E (k)| 

e = "R z max [uTIklt (8.6) 

k=1 


(N = number of control points) . - 

One can see the continuous improvement of the designed wall 
contours. The convergence range decreases continuously during the 
iteration. This phenomenon is to be attributed to the existence of 
harmonics in the approximating u^ distribution. As already dis- 
cussed in Chapter 5.2.1, these contributions can be found by 
iteration using an adapted design rule. However, in practice it 
has been found that already after three to four iteration steps, 
the designed wall shape is sufficient for calculating an adaptive 
contour (see Chapter 8.3 and 9). Therefore, no modification to the 
design rule (5.11) was necessary. 
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8 . 3 Result Comparison with DFVLR Adaptation Program 

As already mentioned in 6. 1.2, the DFVLR has their own 
single step method based on the Fourier method for regulation of 
its three-dimensional adaptive test section (DAM) . It has been 
used successfully. The wall contours calculated with this program 
are therefore "experimentally verified" and can be used to test 
the present adaptation algorithm. 

The two flow cases post calculated with AD API from the "DAM" 
are shown in Figures 15a and 15b together with the DFVLR reference 
values [7, 31]. We show the results for the first adaptation 
step, that is, the wall deflections and pressures calculated based 
on wall pressure distributions measured for aerodynamically straight 
walls „ 


Figure 15a shows the comparison for the FFA model (parabolic 
spindle) for 0° angle of attack and the nominal incident Mach num- 
bers M^ = 0.5. Since this is a flow with rotational symmetry, we 
only show the wall contour and the pressure distribution averaged 
in the circumferential direction. As can be seen, the runs 
calculated with ADAPl agree well with the corresponding reference data. 
The remaining small deviations can be attributed to the additional 
physical correction of the input pressure measured values in the TUB 
method. (see Chapter 7). The effective incident flow Mach number in 


this case is found to be M kor = 0.5036 (blockage 3,1%), 
where the upstream of the model has been removed. 


a value 
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Figure 15b shows corresponding comparison calculations for a 
wing- fuselage combination. This is the AGARD calibration model 
B with 4° incidence angle and M^ = 0.5. The reference values again 
are well reproduced.. The Mach number correction is now AM^ = 0.006 
for a blockage ratio of 3.5%. 
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The good agreement of the calculated results in both flow 
cases therefore clearly verifies the numerical methods of wall 
adaptation developed here, (design algorithm, single step formula) . 
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9 . EXPERIMENTAL RESULTS 
9.1 Tested Models 


Three models were made for testing the 3-D adaptive test section. 

A rotationally symmetric body, the ONERA C5 calibration model, 
is used for pressure measurement. It has 21 static pressure caps 
distributed over the entire model length of 166.27 mm. Figure 16a. 

Its maximum diameter of 24 mm leads to a blockage ratio of 2% in the 
octagon test section. With this model, the first successful 3D 
wall adaptation was achieved at the TUB. [30]. 

The ZKP-F4 model. Figure 16b, has an Airbus-like cross section 
with a transsonic wing (Span 120 mm) . It is designed for 
force measurements which are performed with a three component 
balance installed in the fuselage. Its blockage in the 3D test 
section is 1.2%. 

Figure 16c shows the Canard configuration used for pressure and 
force measurements. The Delta shaped wing has an area which is 
about 4 times larger than the F4 wing with the same span of 120 mm 
(aspect ratio, 2.29) . In this way one obtains not only much larger 
normal forces, but also remarkable risk differences between 
measurements for straight and adapted channel walls. Also large 
changes in the moment variation produced by the effect of the 
Canard control surface, so that the model is especially well suited 
as an indicator for wall interferences. 

The static pressure can be measured from 10 pressure caps 
distributed over the fuselage. The blockage of the model is 1.3% 
in the octagon test section. 
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9.2 Calibration of the 3D Test Section 


Before installing the model, the test section has to be cali- 
brated, that is, the perturbation of the parallel basic flow due 
to the wall boundary layer has to be compensated for by a 
preliminary deflection of the flexible walls. In the case of 
the octagon section, an additional expansion was necessary in 
order to equalize the displacement effect of the model sphere. 

The resulting (aerodynamically straight) wall contour is 
shown in Figure 17 for an incident Mach number of M =0.8. Here 

CO 

we are dealing with a photographic picture of the color monitor, 
which is used for the graphic recording of adjusted wall shapes 
during the tests. The deflections are given in mm. (wall 1-4 top, 
wall 5-8 lower part of the image) . The values are interpolated 
linearly at the actuators over the test section length. 

By means of the same preliminary adjustment shown in the 
circumferential direction, (the deviation in the sphere region is- 
produced by the missing ninth actuator term along the top and 
lower wall) we obtain an almost constant pressure distribution 
at least in the region of the model. Large differences occur 
at the test section beginning before the first actuator because 
of residual wall boundary layer influences. They also occur at /53 

the rear part by the not exactly compensated for displacement 
influence of the sphere. 


In order to assure a reliable wall adaptation later on, we 
decided to only use the wall pressures measured between the first 
and eighth actuator for calculating the external field. The other 
measured values are in place by an extrapolation according to 
Equation 7.10: 


Uj 1 (x) = 


L U i‘ (x D,20) 


X < X 


X > X 


M 


D , 20 


(9.1) 
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The also described correction of the pressure distribution 
for aero dynamically straight test section walls discussed in 
Section 7 did not seem to make sense considering the poor pressure 
measured values at the two test section ends discussed above. 

9.3 Measurement Results with the Canard Model 


In this chapter we will use the example of selected results 
for the Canard model to demonstrate the high performance capability 
of the adaptive octagon test section. An extensive documentation 
of test results is obtained in [32] . 

For the test execution, one important feature was the definition 
of a practical truncation criterion for wall adaptation. One possible 
indicator for this is the convergence behavior of values measured 
on the model. (forces, pressures) or axis - the change in the adjusted 
wall shape in following regulation steps. Of course, eventually /54 
constant wall contour also then leads to constant model measured 
values. However, in the last analysis, removal of the inter- 
ferences is only required in the region of the configuration in the 
flow. Since the wall-induced perturbations towards the mode de- 
creased (M^ < 1) , this has a direct influence on the accuracy limit 
with which the adapted wall shape has to be determined, (see 
Equation 3.4)). This effect is discussed in the next section for the 
special case of a rotationally symmetric flow. The estimation 
obtained there is too "sharp" in practice, because a possible 
extinction of wall interferences in the model region is not 
considered. a basis for evaluating the adjusted wall deflection, 

we therefore use the behavior of the model measured values. 

The typical sequence of wall adaptation is illustrated in 
the following table. 
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Flat 

walls 

1. WK 

2. WK 

3. WK 

C A 

.207 

.159 

.172 

.168 

C W 

.030 

.025 

.027 

.025 

C M25 

.031 

.035 

.034 

.034 


Canard- Model 

M = 0.7, a = 4.54° 

00 


This shows the change in the forces measured on the Canard Model 
at the individual adaptation steps for the incident flow conditions 
M =0.7, a “4.54°. The difference in the coefficients from the 
second to the third wall contour already is within the range of 
the reading accuracy of the balance reading, so that the test section 
walls can be considered to be adapted. 

As a check, we compared the measured wall pressure distribution /55 
after the last regulation step with those calculated for the outer 
field. As discussed in Chapter 3, they have to agree for the 
adapted wall contour. For the flow case under consideration, a 
comparison like this is given in Figure 18 for the upper flexible 
wall. For clarity of the wall interferences, we also show the 
pressure distribution measured for aerodynamically straight walls. 

It can be seen, that the third wall contour canalso be considered 
to be almost adapted even in the "strict" sense. The remaining 
pressure differences in the front area of the model suspension 
apparently are truncation effects, which are caused by the effective 
shortening of the test section along the section between the 
first and eighth actuator (see Chapter 9.2). For the two rear actuators 
9 and 10, we also did not calculate any deflections, but instead 
retained the constant preliminary adjustment used for the prevailing 
incident Mach number . 

As an example for the wall deformations to be adjusted, 

Figure 19 shows the adapted contour for the Canard model at 
M ot = 0.8 and a= 8.04°. The dominating influence of the downwind 
field is especially remarkable, which requires a relatively large 
dropping of the rear test section cross section. (eighth actuator) . 
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The necessary moving back of the wall shape in the forward region 
does not have any important effects on the data measured on the 
model, as the corresponding wind tunnel tests show. 

The next experimental proof of the high quality of the wall 
adaptation is shown in Figure 20. It shows the pitch moment 
polar measured at M m = 0.7 for the Canard model for straight and 
adapted test section roles. For comparison, we show also the 
measured polars in the adaptive profile T2 of the ONERA and the 
polars obtained at the "DAM" for the model. /56 

These measured values can be considered to be free of inter- 
ference because of the very small blockage ratios (DAM 0.06 %, T2 
0.18%) . The T2 tunnel was adapted according to the method developed 
by Wedemeyer (see Chapter 10) . The very accurate reproduction 
of the interference - free moment variation in the octagon test 
section is especially remarkable, because the polar indicates 
strong wall interferences for straight walls. A corresponding 
comparison for the C A (a) and C A (C W ) polar is given in [32]. 

In order to determine the adaptive wall shape, the panel method 
was used for all tests performed - (Program "ADAPl" ) . The average 
calculation time per adaptation step was 4 5 seconds in the wind 
tunnel computer HP 1000F. In order to adjust the test section walls, 
in general 1-2 intermediate steps were required. The adaptation 
was always started with aerodynamically straight walls in order to 
demonstrate the wall interference. It was only at the incident 
Mach number M ot = 0.95 that oscillating wall contours occurred at 
higher incidence angles of the Canard model, so that a reliable 
adaptation could not be reached. This unfavorable regulation behavior 
was certainly attributable to the very poor flow quality in the test 
section for flow near the speed of sound (turbulent separation in the 
nozzle) . 

/ 57 
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On the other hand, the measured wall pressure distribution 
showed local supersonic zones for aerodynamically straight walls. 
This means that the contour determined with the panel method could 
have been wrong. An additional calculation with the TSP adaptation 
program indeed showed somewhat different deflections as is shown 
in Figure 21 for the upper and lower test section wall. 

After improving the flow quality in the test section 
we wish to examine with experiments to what extent these differences 
affect the model flow. 

9 . 4 Influence of the Wall Contour Inaccuracies on the Floor 

on the Model 


In this chapter we will estimate how accurately the test section 
walls have to be adjusted in order to reduce remaining interferences 
in the model measured values through a specified amount. For 
clarity, we will investigate only the effects of perturbations 
in the adapted wall contour which has rotational symmetry 


6 h (x,R) = z A n sin(to n x) , 


_ M n 

“n " n r» 


n = 1,2,.. 


(9.2) 


The perturbation flow induced inside the test section then satisfies 
the two-dimensional potential equation 


+ 7 *r * *rr ’ 0 ’ * 1 - > 0 


(9.3) 


with the boundary condition 

^ r (x,R) = <5h'(x) = X A n U) n COS(c*> n x) 


The solution is 
4> ( x , r ) 


(see Chapter 

= j £ A n cos(u> 


5.2.1) 

x) v-v*-) 

" I i (“ n BR) 


r « R 


(91.4) 


(9.5) 
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For the perturbation speed along the tunnel axis r * 0 
and therefore in the region of the model in the flow we have 


u(x,o) = <fr x (x,o) = - i 2 


„ - sin (a) x ) 0.6) 

6 ii^eR) 


The comparison (9.6) and (9.2) shows that the wall contour 
perturbations decay towards the center of the tunnel, and 
the damping factor is given by 

D(»J = A" 

" e i 1 (« n eR) (9-7) 

X 

Because of ^j(x) > for x > 0 /13/, we find 


D ( in ) i 


(9.8) 


and from this 


U(X,0)U 7D(»„) I A J 4 4 B 2 r ^ 


ElA.I « -?r 'ii E V 0-9) 


where AS is the minimum path transducer preparation. 


Use of the Parseval equation /33/ 


l A 2 $ r f 6 h 2 (x)dx ^ 2 C : 

n l 

0 


for 


(9.10) 


6h(x) | $ C 


(9.11) 


finally gives the simple estimation 


u(x,°) I « — 

6 2 R 


'f c 


(9.12) 
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If the maximum permissible remaining interference in the model 
region is given by 

| u ( X ,0) I * 6U (9.13) 

then from (9.11) and (9.12) we obtain a tolerable contour error 
of 

(i-H» 2 ) R 

< — 

2 

It can be seen that the adapted wall shape has to be adjusted 
more precisely with increasing incident flow mach number, which 
is physically understandable. 

Quantitatively, for example, for M„ = 0.7, 4,u = 0.006 (R = 83 mm, 
AS = 35 mtn, L * 400 mm) we have the accuracy requirement 

l 6 h ( x ) | s 5 * 6u = 0.03 mm (9.15) 

which at least has a realistic order of magnitude. 


I« h ( x ) I 


/ 


AS 

7l 


<5 U 


(9.14) 


10. A NEW METHOD FOR TESTING 3-D MODELS IN WIND TUNNELS 
WITH TWO FLEXIBLE WALLS 

10.1 Introduction 


Wind tunnels with two adaptive walls have been used for about 
ten years successfully for profile tests /4/. Therefore it is 
natural to also use these test facilities for testing three- 
dimensional models. Even though in this case one cannot achieve 
a completely interference free flow field in the test section 
by deflecting the upper and lower tunnel wall (for example, 
matching of porosity) one can attempt to at least reduce wall 
interferences in the model region. 
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This agrees with the idea that in order to avoid blocking 
effects it should be sufficient to produce a stream tube with 
the same cross-sectional area variation as is used for a complete 
three-dimensional wall adaptation. As Wedemeyer demonstrated with 
a numerical simulation, the resulting symmetric wall shape removes 
exactly the blocking speed oh the tunnel axis /24/. 

Based on this result, he developed a single step adaptation 
method. With it, by super position of one deflection and an 
opposite deflection of the two flexible test section walls, 
one can bring about an interference-free flow along the tunnel 
access (x, y = -Jl, 0)# Figure 22. (Symmetric models without slip 
angle are assumed.) The wall induced perturbation speeds u^(x, 0) 

blocking and Vy(x, -jr* (downwind) to be compensated for are 
then calculated using a singularity representation of the channel 
flow from the previously adjusted upper and lower wall contour 
and the prevailing pressure distributions there /25/. 

The first experimental testing of this adaptation method 
was performed in the adaptive profile tunnels of the TUB and 
the ONERA/CERT /4/. A body of rotation at 0° incidence angle and 
a Canard model which produces lift (see Chapter 9) were tested. 

The pressure distribution measured along the top side of the body 
of revolution also agreed very well at high subsonic flow Mach 
numbers (M^ = 0.84) with the corresponding interference-free 
values. This then confirms the possibility of already being able 
to essentially remove the blocking effect with a two-dimensional 
wall adaptation. In the more general case of flows with lift 
it was also possible to achieve a reduction of the perturbation 
speed, see Figure 20. 

Based on the very low blocking ratios of 0.18% assumed, 
these results, however, only give conditional conclusions about 
the quality of the wall deflections calculated with the Wedemeyer 
program. 
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Independent of this, a theoretical analysis of the adaptation 
method in general leads to the conclusion of an increase in the 
remaining interferences with increasing model incidence angle. 

This is primarily due to the fact, that the wall influence is 
removed independent of the actual position of the investigated 
model, always along the tunnel axis. 

In order to overcome this advantage, we will now introduce 
a new adaptation method based on the 3-D single step method, 
with which it is possible to do the following 

(1) to exactly determine the wall influence in the test 
section; 

(2) with which one can achieve an interference-free flow 
along the model axis. 

The calculation of the interferences and the required wall 
displacement will be described in the following sections. 

10.2 Determination of Wall Interferences in the 2-D Test Section 


The wall induced perturbation speeds are obtained immediately 
by using the general representation formula (6.15) for the wall 
dependent potential, derived within the framework of the single 
step method 




1 

4 TV 3 


/ (v A (q) ■ g(p,q)dq 


( 10 . 1 ) 


9 ( P » q ) = tp-(x p -x q ) 2 + (y p -y q ) 2 + ( V 2 q )2r * 


B 2 = 1 -M 2 
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where the control surface again corresponds to the test 
section geometry# Figure 22. 

The normal speed distribution v^. is then determined in a 
known way by the variation of the upper and lower flexible walls 
and the two fixed side walls (Vj = 0) . The V^-component is deter 
mined within the framework of an outer field calculation based 
on the measured wall pressure distribution (see Chapter 5). 

It is assumed that a measurement of the static pressure along 
the central lines of the two flexible walls and along the 
side wall is sufficient. 


By differentiation of the relationship (10.1) we can then 
directly determine the perturbation speeds in the test section 


d <j> 

u w (p) = m (p) 

3 4 > 

v w<P> ■ <P> 

3 d) 

w w (p) . ^ (P) 


( 10 . 2 ) 


without any assumptions about the model-induced flow field 
(Chapter 6.1.2) directly from the measured wall pressures and 
the previously adjusted wall contour. 


This procedure has the following advantages compared with 
the Wedemeyer procedure discussed above: 


a calculation of the channel inner flow is not required. 
In this way, we avoid assumptions about the model-induced 
perturbations. The division into model effect in a 
symmetric displacement and asymmetric lift effect for 
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high model incidence angles used in this method of 
Wedemeyer is certainly only approximately valid. 

The case of already deflected test section walls does not 
require any special calculations. In the Wedemeyer method 
on the other hand, the measured wall pressure distribution 
always has to be "corrected" on (aerodynamically) straight walls. 

10.3 Calculation of the Adapted Wall Shape 

10.3.1 Fundamental Remarks 

By the deformation of the upper and lower flexible walls, 
the three-dimensional model flow in the test section is super- 
imposed with a two-dimensional flow field. Let us assume 
small (additional) deflections 


h 0 (x) - / vU,h)d5 (u ?e er wall > 

(10.3) 

h u (x) = / vU,0)dC (lower wall) 


Then the induced velocity field again can be assumed to be 
irrotational. The corresponding potential <}, ( x ,y ) then satisfies 
the following linear differential equation for subsonic flow. 


3 2 4> xx + *yy = 0 » 6 2 = 1 “ M co 2 (10.4) 

with the boundary conditions defined by (10.3). 

fy (x,h) = v(x,h) = : v Q (x) 

3d, , (10.5) 

ay(x,O)=v(x,0)=:v u (x) 
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The wall displacements h Q (x) and h u ( X ) i win then have to be 
determined so that the speeds which result from (10.4,5) 


u(x,y) = (x,y) 

( 10 . 6 ) 

v(x,y) = ( x ,y ) 

will compensate for the wall interferences (10.2) as much 
as possible in the model region 

i 

= mi n 

(10.7) 

i 

= min 

The components in the span direction w w can apparently 
not be directly reduced. In the following we will therefore 
assume symmetric models without slip angle. 

For the practical application of the 2-D wall adaptation, 
the physical requirements (10.7) however, have to be formulated 
so that the specified values of u and v will lead to an explicit 
calculation of the wall deflections h Q < x ) and h u ( x ) . The 
condition which makes sense for this case is the following 

u w (x,y M ,0) + u(x,y M ) = 0 
v w (x,y M ,0) + v(x,y M ) = 0 

where 

y M (x) = - x tana + -jj 


( 10 . 8 ) 


(10.9) 


( U w + u) I , 
w Model 


(v + v) | 
w Model 
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is the model axis (Figure 22) . 


In the next section we will first calculate the adapted 
wall shape for the special case a = 0 (model with no incidence) . 


10.3.2 Model at 0° Incidence Angle 


For a = 0 we wish to exactly cancel the interferences along 
the tunnel axis by deflecting the two-flexible test section walls 
according to (10.8) 


u(x,£) *-u w (x,£,o) 
v(x,-j) =-V w ( x »Y»°) 


(Blockage) 


(Downwind) 


( 10 . 10 ) 


In order to calculate the required displacement paths 
h Q £x) , h u (x) we will start with the boundary value problem 
(10.4,5) and wish to establish an analytical relationship between 
u(x,£), v(x,J) and v Q (x), v u (x) . 

For this purpose we will transform (10.4,5) using the usual 
coordinate transformation 

x = £ 

« ( 10 . 11 ) 
y = ph 

into the equivalent incompressible flow problem 
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By introducing the stream function ^ 



U 

( 5 ,n) = || 

(e.n) - f| (S.n) 

follows that 
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(10.14) 
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The Dirichlet problem (10.14) occurs in a similar form also in 
connection with two-dimensional wall interference correction 
procedures, so that we can derive the solution from the appropriate 
references /26,27/. For the speeds induced along the tunnel axis 
we obtain the following clear relationships /26/ 
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v u (t) + ''o (t) dt 


2 cosh 
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Pi 


A reverse transformation into the x,y-plane, because of 

u (5, -ft ) = U (X, |) 

(10.16) 

v (5, |t ) = i »(x, | ) 

gives the following coupled integral equations for calculating 
the adaptive wall shape 
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u(x.-jj) = - 
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V 0 (t) ' V u (t) Ht ! 


u w (x ’2’ o) 


(10.17) 


v(x -?> ■ 1 


v o (t)+v u (t) 

2 coshli^ii 


dt 


= - v w (x, ? ,o) 


The numerical solution of (10.17) poses no difficulties and then 
immediately leads to the desired deflections by subsequent integra- 
tion of v and v 
u o 

As expected, in order to compensate for the blocking speed, 
one requires a wall deformation which is symmetric with respect to 
the tunnel action (cross-section change) . In the case of models 
with lift, in addition one has to displace the test section cross- 
section parallel in the rear region (asymmetric deflection) . 

As already mentioned in Section 10.1, this result was already 
achieved before by Wedemeyer /25/. However, his calculation of 
the adapted wall contour in the form of a series expansion is 
less favorable than the exact solution (10.17). 

10.3.3 Models with Incidence 

With increasing incidence angle a the models move away 
from the tunnel central line. Therefore it is more accurate 
to remove the wall interferences along the model axis (see 
(10.8.9) ). The calculation of the deflections required for 
this can then be done in many ways. 

With boundary value problem (10.14) in general, a relation- 
ship is established between a wall contour change and the subsequent 
induced velocities in the test section. As Mokry shows, it can 
be solved in closed form by separating the variables (Fourier 
method) /27/ 

^ ( £ » n ) = f ( £ » ti » , a 2 » •••» b ^ » ^2* •*•) no 18') 
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where and b^ are the Fourier coefficients from the (formal) 
expansions for h Q U) and h u ^) . In order to determine these 
constant desired wall displacements, apparently one can again 
prescribe — and — along the symmetry axis or any other line 

» s a X] 

in the tunnel. This procedure for example has been investigated 
within the framework of a study paper /28/. 

On the other hand , one can start with the exact relationships 
(10.17) for a = 0 and from this obtain an approximate solution 
for a f 0 using a Taylor series expansion for u(x,y^(x)); and 
v(x,y^(x)) • 


u(x,y M (x)) = u(x,£) - x tan a u y (x,£ ) 


+ x 2 1 a nfa u 
2 yy 


(X, Tg) £ ... 


L 

= u(x,-j) - x tan a v x (x,-^ ) 


(10.19) 


$ 2 x 2 tan*g 


xx 


( x , -IjO 1 • • • 


v(x,y M (x)) = v(x,!}) + 6 2 x tan a u x (x >? ) 


8 2 x 2 tan z a 
2 


xx 


(x, -jr) - 


and in order to replace the derivatives with respect to y one uses 
the equations 


S 2 u + v = 0 
x y 


u - v v = 0 

y x 


(continuity) 

(irrotationality) 


( 10 . 20 ) 
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If we now introduce the abbreviations 

v 1 ( t ) = v u (t) - v Q (t) 

v 2 (t) - v u (t) + v 0 (t) 

1 1 

( t-X ) = jTfi - 2 ir ( t -~xT (10.21) 

1 pi 

K 2 (t ' x)= ^ 2 cosi 


then the adapted wall shape 


x 

h 0 ( x ) = \ S { v 2 ( t ) - v^t)} dt 


(10.22) 

X 

h u (x) = \ f (v 2 (t) + v^t)} dt 


finally results by replacing (10.17) in (10.19) 

— + “ d 

/ Vj (t) ' Kj ('t-x)dt - x tan a / v 2 (t) Hx K 2^ t_x ^ dt 


$ 2 x 2 tan 2 a 
2 


+ oo 


/ v^t) 


-gp: K. 1 ( t-x)dt + 


(10.23) 


* - u w (x»y M (x) s o) 

+ °° 

/ V 2 (t) K 2 (t-x)dt + e 2 x tan 

_ ob 


+ oo 

a / 


oo 


v i ( t ( t-x ) d t 


B 2 x 2 tan 2 or 


+ 0° 


f v 2 ( K 2 (t-x)dt 


v w (x»y M (x) s o) 
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For a = 0 (10.23) of course reduces to solution (10.17). 

The number of terms to be considered for a 0 has to be determined 
by experiment. The linear equation system which approximates 
(10.23) always has the form 

A 5 t * B 5 2 = - “w 

; (10.24) 

C * D ; 2 ■ -K 

and can easily be solved by inversion 


1 i: ' 


1 A 


-1 


A CM 
1 > 


i c 



1 W 


(10.25) 


11. SUMMARIES 


The ever increasing requirements for quality of wind tunnel 
results in high subsonic flows led to the development of so-called 
"adaptive" test sections around 1970. One uses the natural idea 
of removing the cause of wall interferences by matching the test 
section walls through the side ways unlimited model flow. 

After the use of this new technology for profile tests, 
in 1979 we started building three-dimensional adaptive test sections. 
Out of the three realized designs, the special construction features 
of the octagon test section with eight flexible walls has been 
discussed in detail which was built at the TU Berlin. 

After this, we gave a detailed discussion of the newly 
developed numerical methods for 3-D wall adaptation developed here. 
Starting with conventional iteration adaptation methods, we first 
investigate the three-dimensional design path (outer field calcula- 
tion) which has to be solved at each regulation step. Using a 
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simple perturbation calculation, the ill-posedness of this 
inverse problem is demonstrated. Then we introduced the design 
method developed to solve this problem. In this, the desired 
wall shape is calculated directly from the prescribed pressure 
distribution, but instead by stepwise modification of an initial 
contour. The required geometric changes are determined using 
a design rule from the difference between the actual and nominal 
pressure distribution. Each of the actual values can be determined 
with the panel or the TSP method as desired. The post calculation 
task can be solved numerically without problems and in the linear 
case can be reduced to a simple matrix vector multiplication. 

The design method therefore does not only allow a simple check 
of the calculated wall deflections, but has addition numerical 
advantages. 

Use of the iteration wall adaptation method with the prescribed 
outer field calculation method in 1982 gave the basic functional 
proof of the octagon test section. But in order to regulate the 
eight flexible walls, in general more than ten intermediate steps 
were necessary even after matching up the regulation factor which 
determines conversion. 

In order to make a systematic investigation of the adaptation 
behavior, we performed a numerical simulation. The wall pressure 
distribution in the individual regulation steps is not determined 
by a wind tunnel test, but using physically reasonable assumptions 
about the given wall contour. This numerical study did not only 
verify the experimental results, but also shows that even if one 
uses an optimum regulation factor, one can not achieve fast 
conversions. 

In a further theoretical analysis, the allowable regulation 
range finally was determined exactly. It is found that the 
iteration 3-D wall adaptation method diverges ;theoretically for 
all regulation factors. The conversions found in practical 


wind tunnel tests or in numerical simulations are only caused by 
the (inaccurate) numerical approximation of the boundary 
conditions I 

On the otherhand, using the relationships derived during 
the simulation for the channel inner flow, we were able to 
establish a basic formula for the wall dependent perturbation 
potential in the subsonic case. This allows a calculation of 
the interferences in the test section based on the known wall 
contour and the prevailing pressure distribution. A representa- 
tion of the model flow is not required. 

Use of this formula leads directly to a single step 
adaptation method. The calculation of the adapted wall contour 
from the wind tunnel measured values is then possibly based on 
a simple .linear relationship. As a comparison with the iteration 
.method shows, the single step formula is found formally by 
replacing the constant regulation factor by a matrix which 

considers the incident flow mach number and the test section 

geometry. In the case of flow near the speed of sound or 
very large wall interferences, the assumptions used in this 
formulation, however, only apply approximately, so that in these 
cases several adaptation steps may be necessary. 

An experimental testing of the single step adaptation method 
in the octagon test section then clearing confirms the assumed 
fast conversions. The wall contour in general was already 
adapted after two regulation steps in general. We use the 

change in the data measured on the model as the truncation 

criterion. The high quality of the adjusted wall shape is 
demonstrated as an example for a tested model (Canard configura- 
tion) for a incident Mach number = 0.7 using interference- 
free comparison measurements. 
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The outer field calculation required in any adaptation 
step was always performed using the subsonic panel method 
during the testing, in order to save calculation time. It 
was only in one flow case, where the local supersonic regions 
extended to the test section walls, that one could not achieve 
a reliable adaptation. A comparison calculation with the 
TSP program then resulted in a somewhat different wall 
deflection for the first adaptation step. After improvement 
of the poor flow quality which occurs for incident flow near 
the speed of sound in the model test section, one should 
then perform a complete wall adaptation with the TSP method. 

At the conclusion of this paper we then developed a new 
method for testing three-dimensional models in wind tunnels 
with two flexible walls. By using the already mentioned 
general representation formula for the wall dependent 
perturbation potential ,one can then calculate the interferences 
in the test section without any additional assumptions about 
the model induced flow field. It is shown that by a.i.Suitable 
deflection of the upper and lower flexible walls, one can always 
remove the perturbation speeds along the model axis . The adapted 
wall shape for models with incidence is then determined from a 
series expansion from the deflections for a - 0 . 
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13. 


APPENDIX 


13.1 Calculation of the Influence Coefficients 
Using the abbreviation 


♦ij = ' Tie / 9( p i > q )<J q 


the following is obtained for the influence coefficients 
(Formulas (5.31) and (5.36) ) in the body-fixed (x,y,z)- 
coordinate system . Figure 9 





n i,y ay 


1 J 


+ n. 


Z 9 Z 





The calculation of is done however/ in the panel-related 

system (x 1 > y > z — I in the most simple way 


z* 



+ a 


/ ~r\ " * 

/^(x*-C) 2 +(y-n) 2 +z 12 

-a 
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The reverse transformation into the body-fixed coordinates 
is then given by 


M = 9 <P 
3X 3X ' 


9y 




3 <j> 

w 


+ n 


y.j 


3 4> 

TF 


3 (ft 
9Z- 


n 


9_£_ 
y,j 9y 


+ n 




3 4> 


Depending on the distance between the target point P and 
the induced panel, after this three different formula collections 
are used for the perturbation speeds . The range of validity 
was limited as follows according to a suggestion by Hess and 
Smith /17/ (see Figure above) 



< 2.45 


2.45 < 



4 4 



> 4 


exact formulas 

raultipole expansion up to including 
quadrupole term 

source approximation 


In the first two cases, the already mentioned transformation 
into the panel coordinate system is required. Exact formulas: 


With the abbreviations 


1 

d = {-J- (x'-a) 2 + (y'+b) 2 + z ,2 } ? 

P 

, 1 
e = (gi (x' + a) 2 + (y' + b) 2 + z' 2 } 1 
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f = ijr (x'-a ) 2 + (y'-b ) 2 + ?’ 2 } 1 


g = (p- (x'+a ) 2 + (y'-b ) 2 + z' 2 } 7 


We find 


3 <J> 

3X 1 

(x 1 

»y* » 

3 <|) 

3X 1 

(0, 

,0,0) 

3 <f> 

sy ' 

(x ' 

,y\ 

3 ( p 
3y 1 

(0, 

,0,0) 

3 <f> 

3Z ' 

(x* 

.y' , 


4tt 8 


y'-b+f 


ln y'+b+e, 
1n y T ^Tg } 


( v i w * \ 1 m ^ x '+ a+ eg 1 x 1 + a + ee ^ 

r lx ,y ,Z ) - J- (In TT- . - aVftf - ln V ' -a+ftH } 


x ' -a+Bf x 1 - a + 8 d 


= sign (z * ) (sign (x '+a)‘ 


x ' +a 

(y'+b) 

1 2 ' 

ee 


arc tg 


sign (x'-a) [arc tg -L^-> 

I ■ IS • a \ 


x'-a 


, _ + „ | x ' -a 
arc tg J — r— r 




x ' +a 

(y'-b) 

2 ' 

eg 


(y *+b) 
“TT 


3 


3 $ 


3z 


t (0,0,0) - 


Multipole expansion: 

Since the coordinate origin is at the panel center of 
gravity, no dipole term results. With the separation 


/(^) 2 + y' 2 + Z' 2 ' 

p 


78 


it follows that 


8 

1Tx~ 


(x' ,y 



(x 


, 4ab 


r 1 , 1 3a 2 +b 2 B 2 

i-’rv^TF rT 5 


5 (ax') 2 + (B 2 by') 2 -u 

" n J } 


a<t> 

ay^ 


(x‘ ,y 



{y ' 4ab 


r 1.1 a 2 + 3b 2 B 2 
L •r7 T+ 7I r rT 5 


5 (ax' ) 2 + (e 2 by ' ) 2 -| i 

" 7F" FT 7 J> 



>y 



a 2 + b 2 B 2 

rT 5 


5 

7 F" 


(ax' ) 2 + ( B 2 by 1 ) 2 j } 


Source approximation: 

In the body-f ixed reference system we assume that the 
panel center of gravity has the coordinates (x 0 »yo»z 0 ) and 
the target point is given by (x,y,z). 

Using 


f T 

To = / - 7 TZ 


(x-Xo ) 2 + (y-yo) 2 + (Z- Zo) ; 


we find that 


rl (x > y>z) - 


1 r 4ab Xp-x 

4ttB *• ~ B 2 " 7 3 


H (x.y.x) = 


|| (x.y.z) = 


4 TT B 
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1 [ 4ab ] 


Influence coefficients A-. for M — > 1 ( 6 — >0) 

1 J 00 

For — > 1 the upstream and downstream effects of the induced 

panels vanish. The speeds induced in the circumferential direction 
are 


3 <p 


in 


(y'-b) 2 +z 12 

Ty T + b) 2 +z* T * 


fyr (0 ,y ' , z.' )= sian (z ') {arc tg - arc tg 

z* * 0 

(0,0,0) = -jr 


13.2 Analytical Investigation of the Conversions of the 
Iterative Wall Adaptation Method 

We will start with the iteration scheme (3.3) 


M - V A t u i (n) 3 


V I (n+1) = v i^ + K (v A^ ‘ Ke (0,1) 


n = 1,2,... (i3.i) 


and we will define the increments 
{Ul ( " t1! = u, (n+1) - u, (n) 


6 V 


(n+1) _ (n+1 ) _ „ (n) _ K / v (n) _ (nK 

‘ V I V I ' N ' V A V I ; 


(13.2) 
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The purpose of the following analysis is to establish a re- 
cursion formula for the SVji . 


* K < v A (n) + V A 
= (1-K)6v, (n+1) 

Using relationships (6 . 12 ) - (6 . 14 ) 
v A [6u, (n+1) ] - Ac = A(v A t«u, 


(n+1 ) \ 

I ' 

[*u I <nt,) ] - v, (n) 

* kv a ^**"1 

we then find 
(n+1) ] - 6Vj (n+1) ) 


6 V 


(n+1) } 


- (13.3) 


(13.4) 


and from this 


A ( V A [6 Ul (n+1) ] = - A6 Vl (n+1) 


(13.5) 


(13.3) and (13.5) give the important intermediate result 


A'6v I (n + 2) = { A 1 - K( A+A 1 ) } 6Vj (n+1) = ( A 1 - KE ) 6 v j ( n+ 1 } , 


(13.6) 


n = 1,2,... 


The iteration rule (13.6) apparently can be carried out exactly 
uniquely when the free constant in the solution collection of 
A 1 is suppressed, i.e. we only consider the subspace 

6v I (n+1) eR 0 -!- = Rt Ce’ > C = konst. , ? > = (1,1 1), n = 1,2,... 

This additional assumption amounts to y.^'e R 0 - because of 

M 

(13.2). This only means that out of the first incident number 
of solutions v A [6Uj( n+1 )] "physically correct" solution 
has to be taken (see Chapter 5.1). 
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Since A' of the subspace Ro 1 can be inverted, using the 

abbreviation A'| = A 1 from (13.6) we find the desired recursion 

R - 

relationship ° 


6v T (n+2) = (E-KA'" 1 )6v I (n+1) = ( E-KA ' " 1 ) n 6Vj 


(13.7) 


n 


1 , 2 ,... 


with which one can then investigate the conversion properties 
of (13.1). Because of 

<$ v j ^ n + 2 ^ — >0 < = > ( E - K A 1 ^) n — >0 (0 Matrix) (13.8) 

then according to a well known mathematical theorem /16/, 

(13.1) converges. Then the following estimation applies for 
the spectral radius of the iteration matrix q = E - KA ' ~ ^ 

e(G) = max | A i | < 1 (13.9) 

Ui«N-1 

The Eigen values *■} of G are 

*i = 1 ~ Kp i , i = 1 , . . . ,N- 1 (13.10) 

where u.. = -J— are the Eigen values of A 1 

According to the theorem of Gerschgorin /16/ we have the 
following inequality for the 

|u i " 1 I' 1 , i = 1 N-1 (13.11) 

(This follows from (6.6), (6.8), (6.9) and the fact that the 

are also Eigen values of A). 

From (13.11) we find Re(u^)>0 and therefore the following 
conversion theorem follows. 
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Let U M = ^ e ( u ^) + ilm(u^) be the Eigen value of A # / 
with the smallest real part. Then the iteration procedure (13.1) 
converges exactly for all K with 

0 < K < 2Re(u M ) (13.12) 


Proof : 

From (13.10) it follows that 


X. |2=X i X i =(1-Ky i )(1-Ky i )=1-2KRe(y.) + K 2 |Vi I 


= 1 - 2K 


Re(v.) 


+ K 2 


^7 


= 1 - 


2KRe(v i )-K 2 


and therefore because of 

lx i l 2 < 1< = > 2KRe(v i ) - K 2 > 0 , 1 = 1, ( 13 . 13 ) 

the theorem is immediately proven. 

I 

The matrix a 1 j is a discrete form of the integral operator 
(6.3) (restricted to the partial space of the functions which 
are not constant on S) . One then has to ask how and 

therefore the magnitude of the permissible regulation range 
2Re(v^) depends on the selected paneling. 

From literature we can derive the general fact that 
Fredholm integral equations of the second kind have an infinite 
number of Eigen values, which can at the most accumulate at the 
origin /20/. In our case the Eigen values (6.3) are also restricted 
(estimation (13.11) applies also for the continuous operator). 
Therefore they have to have an accumulation point ' (theorem of 
Bolzano-Weiers trass) . Therefore we have 

Re(v M ) -> 0 for N -> co (13.14) 
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and therefore using (13.12) we find 


K -> 0 for N -> • (13.15) 

The conversion of the iteration method (13.1) deteriorates 
as the accuracy of the inner and outer field calculation increases. 
It even fails completely in the case of infinitely fine paneling! 

The "conversions" found during practical wind tunnel testing 
or during numerical simulation therefore axe only caused by the 
inaccurate numerical approximation of the boundary conditions! 
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Figure 1. Flow Regions for 3D Model Flow. 
Key: 1) control surface S. 
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Figure 2. Iteration Wall Adaptation -Diagram. 

Key: 1) outer field calculation; 2) wall displacement; 

3) new; 4) yes; 5) no. 
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Figure 3. Calspan Pipe: Velocity Measurement System 
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Figure 4. Octagon Test Section with Canard Model. 



Figure 5. External View of the Octagon Test Section. 
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Figure 7. Design Methods. 


& 

Entwurfs- 

zyklus 


Key: 1) input data; 2) iteration starting values; 3) initial 

contour; 4) yes; 5) no; 6) design; 7) post calculation; 
8) design cycle. 
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h(x)~sin (rx) 




Figure 3. Design Rule. 
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Figure 9. Paneling of Control Surface. 


92 













7 


24 

22 

20 


18-1 » 1 1 u# 1 
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1 REGELFAKTOR K 


Figure 10. Convergence of Iteration Wall Adaptation. 
Key: 1) regulation factor K; 2) iterations. 
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Figure 11. Subsonic Programming System. 


Key: 1) definition of surface points; 2) paneling of control 

surface; 3) panel coordinates; 4) determination of influence 
coefficients; 5) control points; 6) measured values; 

7) component if; 8) adapted values; 9) calculation of 
the adapted wall shape. 
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Figure 12. Adaptation Program ADAP 1. 

Key: 1) matrices; 2) post calculation with panel method; 

3) single step formula. 










WANDKONTUR [mm] 


▼ ADAPTIERT 
A EINSCHRITTVERF 3 


Moo =0.7 


r*Ti] 


5 WANDPOSITION 

Figure 13. Analytical Test Case C5-Body. 

Key: 1) pressure; 2) adaptive; 3) single step method; 

4) wall contour; 5) wall position. 


• • 









WANDPOSITION [mm] 


Figure 15a. FFA-Model . 

Key: 1) wall deflection; 2) wall position; 3) pressure coefficients 

M 00 = 0.4994 (M kor = 0.5036) a = 0° i 









2 WANDPOSITION [mm] 


Figure 15b. AGB-Model 3.5%. 

Key: 1) wall contours; 2) wall position; 3) pressure distribution 


M m = 0.4998 (Mk nr = 0.5060) a = 4' 
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Figure 17. Preliminary Adjustment of Walls (M^ = 0.8). 
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Figure 18. Canard Model. 

Key: 1) pressure coefficient; 2) wall position; 3) flat wall; 

4) adapted wall (3rd WK) ; 5) outer field (adapted wall) . 


M m = 0.7 a = 4.54° 
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Figure 19. Adaptive Wall Contour, Canard Model. = 0.8, a = 8.04 
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Figure 20. Canard Model. 
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Figure 21. Canard Model. 

Key: 1) deflection; 2) wall position; 3) lower wall; 

4) upper wall. 

M 0,,= 0.95; a =3.61° 
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Figure 22. 3D-Model in Test Section with Two Flexible Walls. 
Key: 1) control surface. 


106 


POSTSCRIPT 


This work was done during my activity as a scientific 
worker at the Institute for Aerodynamics in Space Flight at 
the Technical University Berlin. 

I would like to thank Professor Ganzer for scientific 
advice and his interest during the writing of the paper. 

I was stimulated by Professor Renner and Barche and 
thank them for carrying out my examination. 

I would like to thank the members of the working group 
of Professor Ganzer for their collaboration and the friendly 
working atmosphere. 


107 


VITAE 


Rainer Rebstock 


born 



1960 - 1964 

1964 - 1972 

May 1982 

After Winter 
Semester 72/73 

May 1975 

October 1979 

November 1979 - 
March 1986 


Lower schools 

Mathematical natural gymnasium 
Abitur 

Study of Mathematics and Physics at 
the TU Berlin 

Preliminary Diploma Examination 
Diploma 

Scientific employee at the Institute for 
Aerodynamics and Space Flight at the TU Berlin 


108 



